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ABSTRACT 

[Fe-Fe]-hydrogenases  are  a  class  of  metalloenzymes  that  catalyze  the  production  of  H2  from  two 
protons  and  two  electrons.  In  this  work,  we  used  density  functional  theory  (DFT)  calculations  to 
analyze  the  mechanism  of  hydrogen  production,  providing  insight  into  the  role  of  the  intermediates  in 
the  catalysis.  We  also  validated  the  exchange-correlation  functional  applied  within  DFT  for  model 
compounds  of  the  active  site  in  [Fe-Fe]-hydrogenase,  enabling  us  a  reliable  application  for 
understanding  previously  established  hydrogen  production  hypotheses,  as  well  as  providing  a  starting 
point  for  a  future  investigation  of  the  effects  of  the  protein  environment  on  the  catalytic  mechanism  of 
[Fe-Fe]-hydrogenases. 

KEYWORDS:  Hydrogenase;  hydrogen;  density  functional  theory 
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INTRODUCTION 


Recently  Prince  el  al.1  assessed  the  potential  efficiency  of  photobiological  hydrogen  production^  as  a 
renewable  fuel,  regarding  the  maximal  energetic  efficiency  when  using  hydro genases.  It  was  noted  that 
additional  improvements  in  hydrogen  generation  at  the  biochemical  level  are  required.2  In  our  first  step 
towards  understanding  the  mechanism,  we  explored  all-iron  hydrogenases.  The  active  site  domains  of 
the  all-iron  hydrogenases  from  the  sulfate-reducing  microorganism  Desulfovibrio  desulfuricans  (DdH) 
and  from  the  anaerobic  soil  microorganism  Clostridium  pasteurianum  (CpI)4, 5  were  found  to  be  very 
similar/  although  it  has  been  inferred  from  the  manner  in  which  the  crystals  were  obtained  that  the  CpI 
structure  represents  an  oxidized  state  and  that  the  DdH  structure  represents  a  reduced  state/  while 
Nicoiet,  et  al.,7  suggested  that  the  DdH  crystal  represents  a  mixture  of  anaerobically  oxidized  and 
reduced  states.  The  crystal  structures  for  active  sites  in  the  [Fe-Fe]-  and  [Ni-Fej-hydrogenases  have 
also  been  shown  to  be  surprisingly  congruent,7’ 8  reflecting  convergent  evolution.*7  However,  although 
[Fe-Fe]-hydrogenases  and  [Ni-Fe] -hydrogenases  readily  catalyze  the  interconversion  between  the 
hydrogen  molecule  and  two  protons  plus  two  electrons,10  limitations,  such  as  the  sensitivity  to  C^,11’ 12 
still  need  to  be  resolved.  Moreover,  an  understanding  of  the  dependence  of  the  protein  environment  on 
the  catalytic  mechanism  is  still  lacking,  as  recently  reviewed  9 
In  order  to  gain  insight  into  the  catalytic  activity  of  [Fe-Fej-hydrogenases,  a  number  of  theoretical 
density  functional  theory  (DFT)  studies  emerged,13'23  focusing  mostly  on  the  Fe2$2  subcluster  in  the 
active  site,  (abbreviated  [2Fc]h),  w'hile  excluding  the  remaining  subcluster  in  the  active  site,  namely  the 
Fe4S4  cubane  (abbreviated  [4Fe-4S]n,  which  is  linked  to  [2Fe]n  through  a  cysteinyl  sulfur  (see  Figure 
I).  Considering  both  theory  and  experiment,  Armstrong9  outlined  the  relationships  between  different 
states  of  [2Fe]n,  addressing  structure,  mechanisms  of  activation,  and  hydrogen  cycling  (see  Figure  3  in 
reference  9)  as  well  as  posing  several  questions,9  In  this  work,  we  have  modeled  the  intermediates  in 
Armstrong’s  outline9  as  a  starting  point  for  elucidation  of  the  mechanism  and  expanded  on  it  by 
including  an  additionally  oxidized  species.  We  endeavor  here  to  confirm  this  outline  and  provide 
insight  into  the  questions  posed.  We  also  validated  the  x-c  functional  to  be  applied  within  DFT  for  a 
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model  system  for  [Fe2]n  and  for  two  other  model  compounds  for  the  CpI  active  site,  enabling  us  a 
reliable  application.  The  structural  changes  we  calculated  for  the  different  functionals  emphasize  the 
importance  of  applying  the  appropriate  level  of  theory,  and  thus  our  work  will  serve  as  the  starting  point 
for  future  investigations  of  the  effects  of  the  protein  environment  on  the  catalytic  mechanism  of  [Fe-Fe]- 
hydrogenases. 

METHODS 

DFT  calculations  at  the  BP8624’  25/LACV3P**26  and  B3LYP25,27’  28/LACV3P**26  levels  were 
performed  with  the  software  package  Jaguar.29  Note  that  the  two  Fe  atoms  in  [2Fe]n  are  termed 
proximal  and  distal  (Fep  and  Fed),  which  refer  to  their  proximity  to  the  [4Fe-4S]  subcluster  (we  also  use 
subscripts  p  and  d  to  indicate  to  which  iron  atom  ligands  are  coordinated).  The  coordination  geometry 
of  [2Fe]H  with  a  carbonyl  (COb)  bridging  the  Fe  atoms  is  that  of  two  edge-bridged  square  pyramids;  the 
CN,  CO,  and  noncysteinyl  sulfur  atoms  form  the  base  and  the  axial  position  is  occupied  by  the  bridging 
carbonyl,  COb.  In  some  unbridged  structures,  COb  rotates  -180°  around  Fed  but  still  occupying  an  axial 
position,  resulting  in  an  inverted  square  pyramid  about  Fej,  although  the  two  square  pyramids  still  share 
an  edge.  We  used  a  dihedral  angle,  <j)cos  t0  characterize  the  orientation  of  COb  in  [2Fe]n;  <j>co  is  defined 
as  X-Fed-Fep-COb,  where  X  is  the  middle  atom  in  either  the  propane  or  dimethyamine  that  bridges  the 
noncysteinyl  sulfur  atoms. 

First,  three  model  compounds  of  the  active  site  were  considered  in  the  calculations  (A,  B,  and  C  in 
Figure  2):  A  is  (p-PDr)Fe2(CO)3(CN)2(CH3S‘)  (5  in  Ref.  17);  B  is  [Fe2]n  based  on  the  CpI  x-ray 
structure4  with  CH^S"  substituted  for  cysteine  and  the  Fe4S4  cubane  removed;  C  is  [Fe2]n  based  on  the 
CpI  x-ray  structure4  coordinated  to  the  Fe4S4  cubane  through  Cys503.  Initial  construction  of  model 
compounds  B,  C,  and  those  in  Armstrong’s  outline9  (model  compounds  1-8)  used  the  coordinates  for 
residues  580,  581,  and  503,  from  the  CpI  crystal  structure,4  i.e,,  for  [Fe2]n>  the  Fe4$4  cubane,  and  the 
cysteine,  respectively.  Hydrogen  atoms  were  added  as  appropriate  and  the  one  of  the  two  lone  oxygen 
atoms  in  residue  580  was  assumed  to  be  a  water  molecule  ligated  to  Fed.  For  B,  C,  and  1--3,  the 
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following  modifications  were  made:  1)  One  carbonyl  group  coordinated  to  Fep  and  one  coordinated  to 
Fed  were  each  changed  to  cyano  groups,  which  is  consistent  with  spectroscopic  evidence30"32  and  with 
the  presence  of  these  groups  ligated  to  the  Fe  atom  in  the  active  site  of  [Ni-Fe]-hydrogenase;33  2)  The 
two  sulfur  atoms  in  [2Fe]n  were  bridged  with  dimethylamine  (-CH2-NH-CH2-).  (This  bridge  essentially 
replaced  the  second  lone  oxygen  atom  in  the  CpI  crystal  structure4  whose  original  assignment  was  made 
in  lieu  of  unresolved  atoms  bridging  the  sulfur  atoms.)  This  bridging  moiety  has  been  previously 
suggested^4  because  it  can  better  act  as  a  base  than  propane,  which  was  originally  presumed.  It  is  also 
probable  that  the  bridging  moiety  is  the  same  for  both  CpI  and  DdH.9  For  A,  however,  we  used  propane 
as  the  bridging  moiety  so  that  we  could  compare  our  results  to  previous  work.17  Note  that  we  denote  the 
noncysteinyl  sulfur  atoms  in  [2Fe]n  with  subscripts  a  or  b  signifying  that  they  lie  either  above  or  below 
the  plane  defined  by  Fep,  Fed,  and  the  middle  heavy  atom  in  the  bridging  moiety. 

For  4-8,  the  following  additional  modifications  were  made:  3)  Following  Zhou,  et  al.,22  a  CHjSH 
group  was  coordinated  to  each  of  the  three  tri-coordinated  Fe  atoms  in  the  Fe4S4  cubane.  The  CH3SH 
groups  might  serve  to  provide  an  approximate  protein  environment  for  the  Fe4S4  cubane  as  well  as  to 
prevent  coordination  between  the  Fe4S4  cubane  and  the  cysteine,  which  is  not  observed  in  the  crystal;  4) 
The  nitrogen  atom  on  the  cysteine  was  capped  by  a  hydrogen  and  the  carbonyl  was  capped  with  an  NH2, 
i.e.,  the  nitrogen  on  the  next  residue  in  the  sequence  was  included  and  capped  with  two  hydrogens.  The 
peptide  bond  between  the  cysteine  and  residue  504  is  thereby  maintained.  Current  model  compounds 
could  therefore  be  used  in  a  quantum  mechanical/molecular  mechanical  mode!  and  results  could  be 
compared  in  a  straightforward  manner.  For  5,  the  bridging  carbonyl  was  rotated  ~180°  and  the  Fe-C-0 
bond  angle  changed  to  180°.  For  6-8,  one  or  two  hydrogen  atoms  were  added  as  shown  in  Figure  2. 

For  A  and  B,  molecular  charges  of  -I,  -2,  and  -3  were  employed  to  simulate  the  oxidized,  semi- 
reduced,  and  reduced  forms.  For  C,  molecular  charges  of  +1,  0,  and  -1  were  employed  because  of  the 
presence  of  the  [4Fe4S]n  subcluster,  which  was  shown  experimentally  to  have  a  +2  charge.1”  Molecular 
charges  for  1  and  2  are  +1  and  0,  respectively,  and  -1  for  3-8.  Spin  states  for  the  semi-reduced  forms 
for  A,  B,  C,  2,  and  6  are  lA\  all  others  are  0.  Geometry  optimizations  were  performed  using  initial 
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coordinates  as  described  above  but,  in  addition,  starting  from  optimal  geometries  for  compounds  lying 
adjacent  in  Armstrong’s  scheme  (Figure  3  in  reference  9).  Geometry  optimizations  were  performed  with 
nitrogen  in  the  dimethylamine  bridging  moiety  oriented  towards  Fep  and  towards  Fed.  For  open-shell 
systems  both  restricted-spin  and  unrestricted-spin  calculations  were  performed.  Geometries  reported 
here  correspond  to  the  lowest  lying  energy  minima  we  found,  except  when  the  conformations  differ 
substantially  from  the  crystal  structure,  i.e.,  those  structures  that  display  an  RMS  greater  than  1  A.  In 
some  cases,  we  found  isoenergetic  optimal  conformations. 

RESULTS  AND  DISCUSSION 

In  previous  calculations  of  model  compounds  of  the  [2Fe]H,  applications  ranging  from  the  gradient 
generalized  approximation  (GGA)  exchange-correlation  (x-c)  functional  within  DFT  e.g.,  Perdew- 
Wang,21,22  BP86,33  to  the  hybrid  B3LYP  x-c  functional18  were  used.  In  a  study  for  a  related  model 
compound,  it  was  attested  recently  that  B3LYP  would  result  in  poorer  agreement  with  experimental 
structures,''1  even  though  that  group  applied  B3LYP  in  earlier  work. Lv  17  The  geometry  of  the  reduced 
form  of  A  (Table  1)  optimized  at  B3LYP/LACV3P**  shows  that  a  carbonyl  group  (COb)  partially 
bridges  the  Fep  and  Fed  atoms  (see  Figure  3),  which  is  consistent  with  previous  work,17  and  with  the 
asymmetrical  bridging  inferred  in  the  crystal  structure  of  (2Fe]n  in  DdH  (Note:  the  initial  assignment  of 
an  oxygen  atom  or  water  molecule  as  the  bridging  ligand  for  [2Fe]n  in  DdH  was  due  to  a  mixture  of 
reduced  and  oxidized  states  leading  to  the  unresolved  carbon  atom.7)  However,  because  it  has  been 
recently  suggested  that  a  non-hybrid  GGA  x-c  functional  may  be  more  appropriate,30  we  also  optimized 
geometries  at  BP86/LACV3P**,  which  effects  distinct  changes;  in  particular,  a  more  tightly  bridged 
structure  results  (see  Figure  3  and  Table  I),  which  is  inconsistent  with  B3LYP-optimized  geometries 
(either  with  ours  or  w'ith  the  results  of  Bruschi,  et  al.i7)  or  with  the  asymmetrical  bridging  found  in  the 
crystal  structure  of  [2Fe]n  in  DdH,6  The  Fep-Fed  distance  we  find  for  the  reduced  state  of  A  is  very 
close  to  the  value  observed  for  CpI,  even  though  this  crystal  structure  is  thought  to  be  an  oxidized  state. 
The  optimal  molecular  geometry  obtained  for  fully  reduced  A  by  Cao,  et  al., 18  who  employed  B3LYP 
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with  a  mixed  basis  set,  displays  no  carbonyl  bridging  Fep  and  Fed*  Interestingly,  the  optimal  geometry 
for  a  compound  similar  to  A  (but  with  a  CH3SH  ligand  instead  of  CH3S’  and,  therefore,  a  -2  charge 
instead  of  -3)  reported  by  Zhou  et  ah,21  who  used  PW92Jb  with  a  basis  set  of  double-zeta  quality,  has 
values  for  Fed-COb  and  Fep-COb  distances  of  1,746  and  2.921  A.  Indeed,  although  Zhou,  et  ah,22  have 
also  extended  their  calculations  to  include  the  [4Fe-4S]n  subcluster,  it  is  difficult  to  assess  the  accuracy 
of  their  results  because  they  were  obtained  with  a  smaller  basis  set  and  the  PW  x-c  functional.  The 
hybrid  x-c  functional  PBE03;  yields  an  optimal  geometry  more  similar  to  the  semi-bridged  B3LYP 
structure  than  to  the  tightly  bridged  BP86  conformation.  The  sensitivity  of  the  x-c  functionals  is 
demonstrated  by  comparison  of  the  partial  atomic  charges  calculated  using  a  Natural  Bond  Order38 
population  analysis  (see  Table  2). 

The  oxidized  and  semrireduced  forms  of  A  do  not  fit  into  Armstrong’s  scheme9  and  are  only  of 
interest  in  terms  of  a  biomimetic  system  because  A  lacks  a  water  ligand,  as  observed  in  the  CpI  crystal 
structure.  For  semi-reduced  A,  we  obtained  bridged  and  unbridged  optimal  conformations,  in 
agreement  with  Bruschi,  et  aL17  For  B3LYP  and  PBEO,  the  differences  in  energy  are  about  4.5 
keal/mol;  for  BP86,  the  difference  is  about  2  kcal/mol.  Geometrical  parameters  are  presented  in  Table  3 
and  compare  reasonably  well  with  the  CpI  crystal  structure.  The  most  notable  difference  is  that  the  Fep- 
Fej  distances  in  the  unbridged  conformations  obtained  with  B3LYP  and  PBEO  are  0,17  A  and  0.9  A 
longer  than  observed  in  the  CpI  crystal  structure.  The  molecular  geometries  we  obtained  with 
restricted -spin  calculations  are  very  similar  to  the  unrestricted-spin  results  in  Table  3  but  the  restricted- 
spin  geometries  are  6,6  and  7.6  kcal/mol  higher  in  energy  than  the  unrestricted-spin  B3LYP-  and  PBE0- 
optimized  structures,  respectively.  For  all  three  functionals,  <S2>  values  for  the  bridged  conformations 
are  about  0,8,  in  good  agreement  with  the  value  expected  for  one  unpaired  electron,  indicating  only 
minor  spin  contamination.  However,  <S2>  values  are  1.6,  for  the  unbridged  conformations  obtained 
with  B3LYP  and  with  PBEO,  which  indicates  additional  spin  contamination  (one  would  expect  <S">  to 
be  3.75  for  additional  unpaired  electrons),  although  for  BP86,  <S2>  for  both  the  bridged  and  unbridged 
structures  is  less  than  0.8. 
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BP86  geometry  optimizations  for  the  oxidized  form  of  A  (Table  4),  when  starting  with  both  bridged 
and  unbridged  geometries,  resulted  only  in  bridged  structures,  while  applying  B3LYP  and  PBEO  yielded 
both  bridged  and  unbridged  structures  with  the  bridged  structure  lower  in  energy  by  0.81  and  4,2 
keal/mol,  respectively,  for  the  two  functionals.  The  unbridged  B3LYP  and  PBEO  structures  have  large 
Fep-Fed  separations  -  3.198  and  3,078  A,  respectively  -  although  the  axial  carbonyl  is  rotated  to  almost 
the  same  degree  as  in  the  bridged  structures,  i.e.,  4>co”173°.  In  contrast,  Bruschi,  et  al.,17  employing 
B3LYP,  found  only  a  bridged  conformation  for  oxidized  A.  It  appears  that  B3LYP/LACV3P**  level  of 
theory  is  able  to  locate  the  two  local  minima  on  the  potential  energy  surface  for  the  oxidized  form  of  A 
corresponding  to  the  bridged  and  unbridged  conformations,  while  B3LYP  with  the  mixed  basis  set  used 
by  Bruschi,  et  ah,17  or  BP86/LACV3P**  are  not.  It  is  also  notable  that  the  Fep-COb  distance  is  shorter 
than  Fed-COb  in  the  bridged  structures  we  obtained  with  B3LYP,  PBEO,  and  BP86,  in  agreement  with 
the  B3LYP-optimized  geometries  obtained  by  Bruschi,  et  al17 
Model  compound  B  with  molecular  charges  of  -1,  -2,  and  -3  and  model  compound  C  with  molecular 
charges  of  +1,  0,  and  -1  are  analogous  to  the  first  two  compounds  in  Armstrong’s  scheme9  and  one 
additionally  oxidized  form,  i.e.,  compounds  1,  2  and  3.  Molecular  geometries  optimized  at 
B3LYP/LACV3P**  were  all  found  to  display  a  carbonyl  (COb)  bridging  Fep  and  Fed,  in  agreement  with 
Armstrong,9  It  appears  that  the  water  ligand  in  B  and  C  precludes  the  unbridged  structures  observed  for 
A.  Geometrical  parameters  agree  with  the  crystal  structures  (see  supplemental  material)  with  two 
exceptions.  For  the  oxidized  form  of  B,  the  separation  between  Fed  and  the  water  ligand  is  0.9  A  shorter 
than  in  the  Cpl  crystal  structure.  However,  the  first  reduction  increases  Fed-water  separation  by  044  A 
and  the  second  reduction  increases  it  by  an  additional  0.3  A,  which  is  consistent  with  the  departure  of 
the  water  ligand  from  the  active  site  upon  reduction.  For  C,  Fed-water  separation  is  longer  than  that  in 
the  Cpl  crystal  structure  by  0,9  A  for  all  three  charge  states.  In  addition,  compound  C  displays  a  close 
contact  between  the  carbonyl  oxygen  in  the  cysteine  and  an  Fe  atom  in  the  [4Fe-4S]n  subcluster.  This 
contact,  w  hich  is  not  observed  in  the  cry  stal  structures,  results  from  deformation  of  dihedral  angles  from 
the  values  observed  in  the  crystal  and  is  prevented  from  occurring  in  model  compounds  1-8  by  the 
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CHjSH  groups  appended  to  [4Fe-4S]n.  The  frontier  molecular  orbitals  for  the  three  charge  states  of  C 
(see  Figure  4)  are  localized  on  both  the  [4Fe-4S]H  and  [2Fe]n  subclusters  as  well  as  on  the  sulfur  bridge. 
In  addition,  the  total  spin  density  for  semi-reduced  C  resides  mostly  on  [4Fe-4S]H  and  only  partially  on 
[2Fe]H  (see  Figure  4),  in  contrast  to  results  reported  by  Zhou,  et  a!.,21  who  found  that  the  total  spin 
density  resides  on  [2Fe]H  for  a  model  compound  that  excludes  [4Fe-4S]H.  It  is  therefore  obvious  that 
the  [4Fe-4S]n  subcluster  must  be  included  in  model  compounds  for  the  [Fe-Fe]-hydrogenase  active  site 
if  further  mechanistic  details  are  to  be  discerned. 

Geometrical  parameters  for  compounds  1-8  are  presented  in  Table  5.  Agreement  between  the  crystal 
structures  and  the  optimal  geometries  is  good  except  that  the  distance  between  Fed  and  the  water  ligand 
in  the  CpI  crystal  structure  is  0.8  A  longer  than  in  1-3,  which  can  be  explained  by  other  forces  in  the 
crystal  acting  on  the  water  molecule.  Our  results  agree  with  those  reported  by  Schwab,  et  al.,39  who 
reported  for  [2Fe]n  an  Fep-Fed  separation  of  2.55  A  and  Fe-S  distances  in  the  range  of  2.33-2.35  A  for  a 
structure  composed  of  [4Fe-4S]n  and  [2Fe]n  (but  with  a  propane  bridging  the  sulfurs  in  [2Fe]n)  and 
optimized  with  BP86  and  a  triple-zeta  quality  basis  set.  Optimal  conformations  for  1-8  are  shown  in 
Figure  5  in  the  same  relation  to  each  other  as  outlined  by  Armstrong.9  However,  we  found  for  several 
species  other  conformations  that  were  nearly  isoenergetic  or  only  slightly  higher  in  energy  than  the  ones 
displayed  in  Figure  5;  for  3,  4  and  5,  conformations  w'ith  the  dimethyl  amine  bridge  flipped  toward  Fed 
are  only  0.2,  4.6,  and  0.2  kcal/mol  higher  in  energy,  respectively.  For  the  most  part.  Figure  5  follows 
Armstrong’s  scheme  but  the  exceptions  are  notable.  The  optimal  conformation  we  found  for  2  (and  that 
also  displayed  a  small  RMS  difference  with  the  CpI  crystal  structure)  shows  a  carbonyl  bridging  Fep  and 
Fed  but  the  Fep-COb  distance  (1.833  A)  is  shorter  than  the  Fea-COb  distance  (2.457  A).  An  optimal 
conformation  16  kcal/mol  higher  in  energy  has  both  Fe-COb  distances  almost  equal,  about  2.0  A.  We 
found  that  the  equilibrium  between  4  and  5  is  not  between  bridged  and  unbridged  conformations  (as 
displayed  in  Figure  3  in  reference  9);  4  is  semi-bridged,  with  a  Fep-COb  distance  of  2.800  A  and  <bco 
equal  to  173°,  while  in  5,  the  COh  group  has  rotated  away  (d>co=l  34°)  resulting  in  a  Fep-CO  distance 
3.750  A.  However,  the  distance  between  Fej  and  the  carbonyl  coordinated  primarily  to  Fep  is  2.654  A 
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and  the  dihedral  corresponding  to  cj>co  (<|>co'  =  NH-Fep-Fe<i-CO)  is  176°,  so  that  5  can  also  be  considered 
to  be  in  a  semi-bridged  conformation  but  with  Fep  and  Fed  bridged  by  a  different  carbonyl  -  the  one 
coordinated  only  to  Fep  in  4.  Mode!  compound  5  is  lower  in  energy  than  4  by  12.4  kcal/mol. 

The  question  of  whether  the  organic  ligand  bridging  the  sulfurs  in  [2Fe]n,  here  dimethylamine,  acts  as 
a  base  in  the  addition  of  a  hydrogen  atom  to  model  compound  5  was  posed  by  Armstrong.9  Although 
we  did  indeed  find  an  optimal  conformation  for  model  compound  6  in  which  the  dimethylamino  group 
is  protonated,  it  is  65  kcai/mol  higher  in  energy  than  the  lo west-energy  conformation  we  found.  The 
energetically  preferred  location  of  the  first  added  hydrogen  atom  in  Armstrong’s  scheme9  is  the  nitrogen 
of  the  cyano  group  coordinated  to  Fed  (see  Figure  5).  We  also  found  for  model  compound  6  an  optimal 
conformation  with  a  hydride  bridging  Fep  and  Fed  but  this,  too,  is  about  65  kcal/mol  higher  in  energy 
than  the  lowest  energy  conformation,  For  model  compound  7,  formed  when  a  second  hydrogen  atom  is 
added,  we  considered  all  combinations  of  two  hydrogen  atoms  attached  to  the  cyano  ligands  (to  both  the 
carbons  and  the  nitrogens),  to  the  nitrogen  in  the  dimethylamino  bridge,  and  to  Fep  and  Fed  via  a  hydride 
■  bridge.  As  shown  in  Figure  5,  the  lowest  energy  conformation,  labeled  7a,  has  a  hydride  bridge  and  a 
hydrogen  atom  bonded  to  the  nitrogen  of  the  cyano  ligand  coordinated  to  Fep.  A  conformation  9.2 
kcal/mol  higher  in  energy  has  a  hydride  bridge  and  a  hydrogen  atom  bonded  to  the  nitrogen  of  the 
dimethylamino  bridge,  although  it  is  unclear  how  this  conformation  is  involved  in  the  formation  of 
molecular  hydrogen.  Another  higher-energy  conformation  that  does  appear  to  be  involved  in  FF 
formation  is  labeled  7b  in  Figure  5;  it  is  similar  to  model  compound  7a  except  that  the  hydrogen 
originally  bonded  to  the  nitrogen  atom  in  the  cyano  ligand  coordinated  to  Fep  has  migrated  to  the  cyano 
ligand’s  carbon  atom  and  is  separated  from  the  hydride  by  only  1.7  A.  Model  compound  7b  is  20.4 
kcal/mol  higher  in  energy  than  7a  but  only  2.3  kcal/mol  lower  in  energy  than  model  compound  8,  which 
has  H?  coordinated  to  Fep.  A  conformation  2.6  kcal/mol  higher  in  energy  than  8  has  the  NH  in  the 
dimethyl  amine  flipped  towards  Fe,j.  Zhou,  et  al.,22  using  the  same  model  compound,  obtained  a 
structure  similar  to  7b  but  only  for  the  oxidized  and  semi-reduced  states,  i.e.,  with  charges  of  +  1  and  0, 
and  did  not  find  an  optimal  geometry  with  IF  coordinated  to  Fep,  as  in  our  model  compound  8,  for  any 
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charge  state.  Superimposing  either  model  compound  7b  or  8  with  the  CpI  crystal  structure  for  [2Fe]n, 
[4Fe-4S]n,  and  the  cysteine  linking  them,  we  found  RMS  differences  of  0.74  A  and  0.66  A, 
respectively.  (For  this  calculation,  we  changed  two  of  the  CO  groups  in  the  CpI  structure  to  CN  and  we 
included  only  atoms  that  were  common  to  both  the  crystal  structure  and  the  model  compounds.  In 
addition,  we  did  not  include  the  bridging  CO  in  the  crystal  structure  and  the  corresponding  CO  in  the 
model  compounds,  which  is  rotated  ~3  80°.).  Superposition  of  8  with  the  CpI  crystal  structure4  is  shown 
in  Figure  6.  The  largest  difference  is  in  the  shape  and  orientation  of  the  [4Fe-4S]n  subcluster,  which  is 
much  closer  to  a  perfect  cube  in  8  than  in  the  crystal  structure.  The  difference  in  orientation  is  due  to 
coordination  between  the  cycsteine’s  carbonyl  group  and  one  of  the  CH3SH  groups. 

Plots  of  the  highest  occupied  and  lowest  unoccupied  molecular  orbitals  for  7b  and  8  (see  Figure  7)  are 
similar  for  the  two  species;  the  HOMO  is  localized  on  the  [4Fe-4S]n  and  the  LUMO  on  the  [2Fe]n. 
Using  NBO  analysis,38  we  found  that  for  model  compound  7b  the  partial  charge  on  the  hydride  bridging 
Fep  and  Fed  is  -0.06  while  the  partial  charge  on  the  hydrogen  bound  to  CNP  is  +0,27.  For  model 
compound  8,  the  partial  charge  on  the  hydrogen  bound  to  CNp  is  +0.21,  which  is  similar  to  mode! 
compound  7b.  and  the  partial  charge  for  the  other  hydrogen  atom  in  Hi  is  +0.04.  We  removed  the  H2 
from  model  compound  8  and  re-optimized  the  geometry,  hoping  to  re-ohtain  model  compound  5.  The 
resulting  geometry  was  very  close  to  model  compound  5  except  for  a  closer  than  usual  contact  (3.734 
A)  between  a  sulfur  in  [2Fe]n  and  the  Fe  atom  in  [4Fe-4S]n  that  is  coordinated  to  the  cysteinyl  sulfur. 
Imposing  our  requirement  that  conformations  must  be  reasonably  close  to  the  crystal  structure,  we 
rejected  this  conformation  but  do  not  believe  that  it  is  evidence  that  model  compounds  8  and  5  do  not 
represent  species  that  are  in  equilibrium. 

CONCLUSIONS 

By  comparing  calculated  to  experimental  structures  for  [2Fe]n,  we  determined  an  appropriate  level  of 
theory  for  investigating  the  mechanism  of  hydrogen  production  in  [Fe-Fe]-hydrogenase.  In  addition,  we 
showed  that  model  compounds  for  the  [Fe-Fe]-hydrogenase  active  site  should  include  the  [4Fe-4S]H 
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cluster  in  addition  to  [2Fe]M,  although  almost  all  previous  theoretical  studies  include  only  [2Fe]H.  We 
also  showed  that  optimal  geometries  for  some  of  the  intermediates  in  a  previous  hypothesis  by 
Armstrong9  are  somewhat  different  than  previously  thought,  while  for  others,  in  particular  fully  reduced 
species  with  transient  hydrogen  species,  ours  is  the  first  theoretical  confirmation.  We  have  established  a 
starting  point  for  a  future  investigation  of  the  effects  of  the  protein  environment  on  the  catalytic 
mechanism  of  [Fe-Fe]-hydrogenases. 
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Figures  Captions; 

Figure  1 .  The  [Fe-Fe]-hydrogenase  active  site  from  the  Cpl  crystal  structure,  which  consists  of  an  Fe2S2 
subcluster,  termed  [2Fe]H?  linked  to  a  [4Fe-4S]n  subcluster  through  a  cysteinyl  sulfur.  Here,  the 
cysteine  is  included.  Red:  oxygen;  yellow:  sulfur;  blue:  nitrogen;  aqua;  iron;  and  grey:  carbon. 

Figure  2.  Model  compounds  A,  B?  and  C 

Figure  3.  Minimum  energy  conformations  for  the  reduced  form  of  A,  Conformations  for  model 
compound  A  optimized  at  B3LYP/LACV3P**  (a)  and  at  BP86/LACV3P**  (b).  Conformation  (b)  is 
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more  tightly  bridged  with  a  Fep-CO  distance  of  2.181  A,  compared  to  2.730  A  for  (a).  Red:  oxygen; 
yellow;  sulfur;  blue;  nitrogen;  aqua:  iron;  grey;  carbon;  and  white:  hydrogen. 

Figure  4.  Frontier  orbitals  and  total  spin  denisty  for  model  compound  C.  The  IIOMO  and  the  LUMO 
are  displayed  in  (a)  and  in  (b),  respectively,  for  the  oxidized  form  (a  molecular  charge  of  +  1)  of  model 
compound  C.  The  HOMO  and  LUMO  for  the  partially  reduced  form  (a  molecular  charge  of  0)  of  C, 
calculated  with  unrestricted-spin,  are  shown  in  (c)  and  (d),  respectively.  The  HOMO  and  LUMO  for 
the  reduced  form  (a  molecular  charge  of  -1)  are  shown  in  (e)  and  (f).  Total  spin  density  is  plotted  for 
partially  reduced  C  in  (g).  Color  legend  for  atoms  is:  red:  oxygen;  yellow:  sulfur;  blue:  nitrogen;  aqua: 
iron;  and  grey:  carbon;  and  white:  hydrogen. 

Figure  5.  Optimal  conformations  of  model  compounds  1-8  inserted  into  Armstrong’s  scheme9 

Figure  6.  Superposition  of  mode!  compound  8  writh  the  II-cluster  from  the  CpJ  crystal  structure. 
RMS=0.740  A.  Only  common  atoms  were  used  and  are  shown,  as  described  in  the  text.  Red:  oxygen; 
yellow:  sulfur;  blue:  nitrogen;  aqua:  iron;  and  grey:  carbon. 

Figure  7.  Frontier  orbitals  for  model  compounds  7b  and  8:  a):  HOMO  for  7b;  b)  LUMO  for  7b;  c) 
HOMO  for  8;  d)  LUMO  for  8.  Color  legend  for  atoms  is:  red:  oxygen;  yellow:  sulfur;  blue:  nitrogen; 
aqua:  iron;  and  grey:  carbon;  and  white:  hydrogen. 


18 


Table  L  Geometrical  Parameters  for  the  Reduced  Form  of  Model  Compound  A,  a  Compound  Similar 
to  A,31  and  the  Crystal  Structure  DdH3 


Geometrical 

Parameter 

B3LYP 

BP86 

PBE0 

Brusehi,  et 
al.17 

Cao,  et 

al.18 

Zhou,  et 
al2' 

DdH3 

Fe-Fe 

2.612  A 

2.594  A 

2.543  A 

2.655  A 

2.655  A 

2.571  A 

2.566  A 

Fep-CNp 

1.956 

1 .903 

1.923 

1.939 

1.941 

1.88 

1.845 

Fep-COp 

1.748 

1.720 

1.727 

1.733 

1.730 

1.74 

1.712 

Fep-SCHj 

2.451 

2.440 

2.382 

2.406 

2.415 

2.19 

2.505 

Fed-CNd 

1.947 

3.899 

1.911 

1.933 

1.934 

1.88 

1.865 

Fej-COj 

1.742 

1.714 

1.717 

1.725 

1.727 

1.74 

1.797 

Fed-COb 

1.758 

1.802 

1.740 

1.754 

1.794 

1.746 

2.494c 

Fep-COb 

2.730 

2.181 

2.559 

2.638 

2.921 

2,807c 

Fep-Sa 

2.384 

2.396 

2.346 

2.387 

- 

2.28 

2.304 

Fep-Sb 

2.402 

2.410 

2.355 

2,396 

2.400 

2.28 

2.290 

Fed-Sa 

2.408 

2.384 

2.353 

2.404 

2.403 

2.28 

2.279 

Fej-Sb 

2.400 

2.346 

2.351 

2.390 

" 

2.28 

2.298 

Fej-C-0 

167° 

149° 

165° 

-164° 

— 

— 

“ 

<t>coa 

-156 

-171 

-159 

-157b 

“ 

— 

t74oc 

a<J>co  defined  in  the  text. 


bAs  specified  by  Brusehi,  et  aL,17  this  value  was  measured  as  the  rotation  of  the  CO  group  with 
'J  respect  to  the  plane  defined  by  Fep,  Fed,  and  the  midpoint  between  Sa  and  Sb,  which  we’ve  found  to  be 
nearly  identical  to  our  definition, 

cThe  oxygen  atom  originally  assigned  as  the  bridging  moiety  is  the  DdH  crystal  structure  was  used  in 
lieu  of  COb-  Note:  the  initial  assignment  of  an  oxygen  atom  or  (water  molecule)  as  the  bridging  ligand 
for  [2Fe]n  in  DdH  was  due  to  a  mixture  of  reduced  and  oxidized  states  leading  to  the  unresolved  carbon 
atom.7 
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Table  2.  Partial  Charges  from  a  Natural  Bond  Order38  Population  Analysis  for  the  Reduced  Form  of  A 
Calculated  with  Different  Functionals 


Bruschi,  et  al.17 

B3LYP 

BP86 

PBE0 

Fep 

-0.022 

0.0045 

-0.276 

-0.126 

Fed 

-0.249 

-0.263 

-0.343 

-0.352 

SCH3 

-0.530 

-0.613 

-0.548 

-0.575 

S2C3Hs 

-0.612 

-0.689 

-0.525 

-0.601 

cob 

-0.208 

-0.145 

-0.154 

-0.137 

cod 

-0.204 

-0.153 

-0.103 

-0.128 

CNd 

-0.565 

-0.565 

-0.536 

-0.549 

COp 

-0.019 

0.0200 

0.0189 

0.0487 

CNp 

-0.591 

-0.601 

-0.535 

-0.582 
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Table  3.  Geometrical  Parameters  for  the  Semi-Reduced  Form  of  A  and  the  Crystal  Structure  CpI4, 


Bridged 

Unbridged 

Geometrical 

Parameter 

B3LYP 

BP86 

PBEO 

Bruschi, 
et  al.17 

B3LYP 

BP86 

PBEO 

Bruschi, 
et  al.17 

CpI4 

Fe-Fe 

2.592 

2.563 

2.532 

2.597  A 

2.791  A 

2.538  A 

2.706  A 

2.587  A 

2.617  A 

Fep-CNp 

1.943 

1.903 

1.911 

1.927 

1.985 

1.924 

1.960 

1.958 

1.857b 

Fep-COp 

1.766 

1.740 

1.744 

1.768 

1.779 

1.741 

1.759 

1.757 

1.754 

Fep-SCHj 

2.389 

2.370 

2.341 

2.337 

2.402 

2.288 

2.355 

2.335 

2.378 

Fed-CNj 

1.946 

1.912 

1.919 

1.916 

1.943 

1.909 

1.912 

1.930 

1 .863 11 

Fed-COd 

1.766 

1.740 

1.744 

1.772 

1.786 

1.754 

1.763 

1.770 

1.831 

F  e^-C  Ob/axi^i 

1.930 

1.926 

1.927 

1,878 

1.793 

1.747 

1.770 

1.750 

2.043 

Fep-COb 

2.082 

1.976 

2.003 

2.116 

— 

— 

2.100 

Fep-Sa 

2.405 

1.363 

2.357 

2.379 

2.418 

2.324 

2.383 

2.372 

2.319 

Fep-Sb 

2.410 

2.388 

2.360 

2.386 

2.354 

2.346 

2.314 

2.355 

2.340 

Fed-Sa 

2.380 

2.348 

2.336 

2.375 

2.361 

2.304 

2,314 

2.312 

2.337 

Fea-Sb 

2.363 

2.321 

2.319 

2.371 

2.376 

2.323 

2.331 

2.355 

2.320 

Fed-Cb-Ob 

142° 

138° 

140° 

— 

— 

138° 

(t>COJ 

172 

173 

173 

-9° 

-13° 

>9° 

— 

178c 

is  defined  in  the  text. 

bThe  CN  ligands  are  assigned  as  CO  ligands  in  the  CpI  crystal  structure. 

cFor  the  CpI  active  site,  the  midpoint  between  the  two  sulfur  atoms  was  used  to  define  c|>co 
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Table  4.  Geometrical  Parameters  for  the  Oxidized  Form  of  A  and  the  Crystal  Structure  Cpl4 


Bridged 

Un  bridged 

Geometrical 

Parameter 

B3LYP 

BP86 

PBE0 

Bruschi, 

etal.17 

B3LYP 

PBE0 

Cpl4 

Fe-Fe 

2.519  A 

2.507  A 

2.468  A 

2.516  A 

3.198  A 

3.078  A 

2.617  A 

Fe'p-CNp 

1.936 

1.902 

1.906 

1.923 

1.932 

1.908 

1.857  s 

Fep-COp 

1.802 

1.769 

1.776 

1.789 

1.768 

1.749 

1.754 

Fep-SCH3 

2.325 

2.344 

2.287 

2.312 

2.297 

2.263 

2.378 

Fed-CNj 

1.927 

i  .881 

1.902 

1.914 

1.927 

1.901 

1 ,863b 

Fed-COd 

1.811 

1.762 

1.786 

1.797 

1.818 

1.790 

1.831 

Fed-COb/ axial 

2.053 

2.056 

2.046 

2.004 

1.774 

1.756 

2.043 

Fep-COb 

1.911 

1.857 

1.866 

1.929 

3.315 

3.193 

2.100 

Fep-Sa 

2.370 

2.353 

2.324 

2,367 

2.355 

2.299 

2.319 

Fep-Sb 

2.371 

2.333 

2.325 

2.352 

2.402 

2.337 

2.340 

Fed-Sa 

2.335 

2.274 

2.298 

2.338 

2.362 

2.322 

2.337 

Fed-Sb 

2.333 

2.213 

2.298 

2.339 

2.384 

2,349 

2.320 

Fej-Cb-Ob 

134° 

128° 

133° 

177° 

177° 

138° 

<|>coa 

174 

172 

175 

173 

173 

178e 

*§co  is  defined  in  the  text, 

bThe  CN  ligands  are  assigned  as  CO  ligands  in  the  Cpl  crystal  structure.  4 

LFor  the  Cpl  active  site,  the  midpoint  between  the  two  sulfur  atoms  was  used  to  define  <J>co- 
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*'  Table  5a.  Geometry  Parameters  for  Compounds  1-6 


Bond 

Cpl 

1 

2 

3 

4 

5 

6 

Fe-Fe 

2.617 

2.592 

2.710 

2,552 

2.756 

2.651 

2.565 

Fep-CNp 

1.857b 

1.933 

1.926 

1.934 

1.933 

1.917 

1.954 

Fep-COp 

1.754 

1.800 

1.806 

1.795 

1.788 

1.766 

1.806 

FCp"Scy5 

2.378 

2.312 

2.423 

2.345 

2.309 

2.191 

2.425 

Fed-CNd 

T8635- 

1.929 

1.954 

1.945 

1,912 

1.940 

1.815 

Fea-COd 

1.831 

1.842 

1.821 

1.817 

1.814 

1.814 

1.800 

Fed'COb/axial 

2.043 

1.869 

2.457 

1,932 

1.774 

1.761 

1.804 

Fed-OH2 

2.959 

2.125 

2.192 

2,112 

- 

Fep-COb 

2.100 

2.273 

1.833 

2.098 

2.800 

3.750 

Fep-Sa 

2.319 

2.332 

2.394 

2.351 

2.335 

2.367 

2.381 

Fep-Sb 

2.340 

2.346 

2.416 

2.343 

2.322 

2.367 

2.346 

Fed-Sa 

2.337 

2.381 

2.335 

2.357 

2.347 

2.281 

2.350 

Fed-Sb 

2.320 

2.365 

2.341 

2.349 

2.360 

2.291 

2.362 

Scys“Fc  cubarte 

2.321 

2.176 

2.372 

2.286 

2.232 

2.273 

2.425 

FCd"Cb-0|tf axial 

138° 

155° 

121° 

146° 

172° 

178° 

178° 

178c 

174 

-176 

172 

173 

6d 

17 

a<t>co  is  defined  in  the  text. 

bThe  CN  ligands  are  assigned  as  CO  ligands  in  the  CpI  crystal  structure. 

cFor  the  Cpl  active  site,4  the  midpoint  between  the  two  sulfur  atoms  was  used  to  define  <^co- 

dThe  dihedral  N-Fep-Fed-COp  is  176°. 
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Table  5b.  Geometry  Parameters  for  Model  Compounds  7,  8,  and  for  Two  Similar  Model  Compounds 
due  to  Zhou,  et  al.22 


Bond 

7a 

7b 

8 

la[l]oxD 

la[0]sB 

Fe-Fe 

2.578 

2.595 

2.672 

2.564 

2.577 

Fep-CNp 

1.820 

2,029 

1.932 

1.926 

1.912 

Fep-COp 

1.800 

1.817 

1.812 

1.768 

1.765 

FCp*Scys 

2.357 

2.387 

2.380 

2.269 

2.283 

Fed-CNd 

1.931 

1.934 

1.935 

1.864 

1.870 

Fed-COd 

1.788 

1.797 

1.813 

1.779 

1.768 

Fed-COaxial 

1.786 

1.785 

1.762 

1.764 

1.756 

Fep-Sa 

2.383 

2.336 

2.387 

2.263 

2.278 

Fep-Sb 

2.396 

2.393 

2.423 

2.284 

2.294 

Fed-Sa 

2,357 

2.338 

2.310 

2.243 

2.249 

Fea-Sb 

2.359 

2.352 

2.329 

2.252 

2.260 

Scys“F^  cubane 

2.325 

2.327 

2.317 

2.158 

2.165 

Fep-H’ 

- 

1.617 

1.574 

1.635 

1.599 

Fep-H” 

1.630 

1.633 

1.687 

1.642 

1.633 

Fed-H” 

1.656 

1,678 

1.902 

1.681 

1.681 

H’-H” 

4,180 

1.699 

0.882 

1.713 

1.685 

NC-H’ 

— 

1.179 

1.874 

1.207 

1.234 

+coa 

11° 

IT 

12° 

10° 

10° 

a(j>co  is  defined  in  the  text. 

^Notation  used  for  the  model  compounds  in  the  last  two  columns  is  the  same  as  in  reference.""  la[l]ox 
and  la[0]s  are  oxidized  and  semi-reduced  and  have  charges  of  +  1  and  0,  respectively. 
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Figure  i.  The  lFe-Fe]-hydrogenase  active  site  from  the  CpI  crystal  structure  ,  which  consists  of  an 
Fe2S2  subclustcr,  termed  [2Fe]n,  linked  to  a  [4Fe-4S]n  subcluster  through  a  cysteinyl  sulfur.  Here,  the 
cysteine  is  included.  Red:  oxygen;  yellow:  sulfur;  blue:  nitrogen;  aqua:  iron;  and  grey:  carbon, 
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Figure  2.  Mode!  Compounds  A.  B,  and  C. 
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Figure  3.  Minimum  energy  conformations  for  the  reduced  form  of  A.  Conformations  for  model 
compound  A  optimized  at  B3LYP/LACV3P**  (a)  and  at  BP86/LACV3P**  (b).  Conformation  (b)  is 
more  tightly  bridged  with  a  Fep-CO  distance  of  2.181  A,  compared  to  2.730  A  for  (a).  Red:  oxygen; 
yellow:  sulfur;  blue:  nitrogen;  aqua:  iron;  grey:  carbon;  and  white:  hydrogen. 
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Figure  4 
a) 
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Figure  4 
e) 


f) 
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Figure  4 
g) 


Figure  4.  Frontier  orbitals  and  total  spin  density  for  model  compound  C,  The  HOMO  and  the  LUMO 
are  displayed  in  (a)  and  in  (b),  respectively,  for  the  oxidized  form  (a  molecular  charge  of  +1)  of  model 
compound  C.  The  HOMO  and  LUMO  for  the  partially  reduced  form  (a  molecular  charge  of  0)  of  C, 
calculated  with  unrestricted-spin,  are  shown  in  (c)  and  (d),  respectively.  The  HOMO  and  LUMO  for 
the  reduced  form  (a  molecular  charge  of  -1)  are  shown  in  (e)  and  (f).  Total  spin  density  is  plotted  for 
partially  reduced  C  in  fg).  Color  legend  for  atoms  is:  red:  oxygen;  yellow:  sulfur;  blue:  nitrogen;  aqua: 
iron;  and  grey:  carbon;  and  white:  hydrogen. 
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Figure  5.  Optimal  conformations  for  1-8  shown  in  the  same  relation  to  each  other  as  outlined  by 
Armstrong  9 


32 


Figure  6.  Superposition  of  mode!  compound  8  with  the  H-cluster  from  the  CpI  crystal  structure.  RMS  = 
0.740  A.  Only  common  atoms  were  used  and  are  shown,  as  described  in  the  text.  Red:  oxygen;  yellow: 
sulfur;  blue:  nitrogen;  aqua:  iron;  and  grey:  carbon. 
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,  Figure  7 

i  a) 
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Figure  7 

c) 


d) 


Figure  7.  Frontier  orbitals  for  model  compounds  7b  and  8:  a):  HOMO  for  7b;  b)  LUMO  for  7b;  c) 
HOMO  for  8;  d)  LUMO  for  8.  Color  legend  for  atoms  is:  red:  oxygen;  yellow:  sulfur;  blue:  nitrogen; 
aqua:  iron;  and  grey:  carbon;  and  white:  hydrogen. 
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Supporting  Information  Available 


Geometry  Parameters  for  model  compounds  B  and  C  are  included  as  Supporting  Information. 
Table  SI:  Geometry  Parameters  for  Model  Compound  B 


Bond 

Cpl4 

Oxidized 

Semi-Reduced 

Reduced 

DdHJ 

Fe-Fe 

2.617  A 

2.558  A 

2.595  A 

2.597  A 

2.566  A 

Fep-CNp 

1.8571’ 

1.929 

1.941 

1.950 

1.845 

Fep-COp 

1.754 

1.798 

1.769 

1.746 

1.712 

Fep-SCFh 

2.378 

2.360 

2.391 

2.454 

2.505 

Fep-CNd 

1.863'1 

1.941 

1.942 

1.941 

1.865 

Fea-COd 

1.831 

1.805 

1.769 

1.736 

1.797 

Fed-COb 

2.043 

2.177 

1.933 

1.763 

2.494“  ~ 

Fed-OH2 

2.959 

2.050 

3.414 

3.741 

Fep-COb 

2.100 

1.883 

2.077 

2.494 

2.807“  " 

Fep-Sa 

2.319 

2.371 

2.404 

2.399 

2.304 

Fep-Sb 

2.340 

2.383 

2.409 

2.417 

2.290 

Fed-Sa 

2,337 

2.372 

2.387 

2.413 

2.279 

Fej-Sb 

2.320 

2.340 

2.356 

2,385 

2.298 

Fed-Cb-Ob 

138° 

131° 

142° 

161° 

“ 

<t>COC 

178 

176 

173 

-162 

1 74d 

a<j>co  is  defined  in  the  text, 

bThe  CN  ligands  are  assigned  as  CO  ligands  in  the  CpI  crystal  structure. 

cFor  the  Cpl  active  site,11  the  midpoint  between  the  two  sulfur  atoms  was  used  to  define  (]>co- 

dThe  oxygen  atom  assigned  as  the  bridging  moiety  is  the  DdH  crystal  structure  was  used  in  lieu  of 
COb. 
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Table  S2:  Geometry  Parameters  for  Model  Compound  C 


Bond 

CpI4 

Oxidized 

Semi-Reduced 

Reduced 

DdHJ 

Fe-Fe 

2,617 

2.634  A 

2.684  A 

2,542  A 

2.566  A 

Fep-CNp 

l.857b 

1.937 

1,942 

1.939 

1.845 

Fep-COp 

1.754 

1.822 

1.817 

1.806 

1.712 

Fep-Scysieinc 

2.378 

2.261 

2.283 

2.342 

2.505 

Fed-CNd 

1.863s 

1.922 

1.927 

1.940 

1,865 

Fed-COd 

1,831 

1.848 

1.832 

1.814 

1.797 

Fed-COb 

2.043 

1.836 

1.841 

2,066 

2.494d 

Fed-OH2 

2.959 

2.042 

2.051 

2.044 

— 

FeP'COb 

2.100 

2.438 

2.459 

1.956 

2.807d 

Fep-Sa 

2.319 

2.317 

2.331 

2.352 

2.304 

Fep-Sb  . 

2.340 

2.322 

2.332 

2.370 

2.290 

Fed-Sa 

2.337 

2,412 

2.421 

2.376 

2.279 

Fej-Sb 

2.320 

2.382 

2.387 

2.346 

2.298 

Fed-Cb-Ob 

138° 

162° 

162° 

137° 

— 

<t>coa 

1 7  8C 

175 

175 

178 

174d 

a<t>co  is  defined  in  the  text. 

bThe  CN  ligands  are  assigned  as  CO  ligands  in  the  CpI  crystal  structure. 

cFor  the  CpI  active  site,4  the  midpoint  between  the  two  sulfur  atoms  was  used  to  define  4>co- 

dThe  oxygen  atom  assigned  as  the  bridging  moiety  is  the  DdH  crystal  structure  was  used  in  lieu  of 
COb. 
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